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ABSTRACT 


Mice oerumetinewOtmbarotropically unstable disturbances in 
mew enOptes was examined with a two-level quasi-geostrophic 
model. The forecast equations were linearized and non- 
dimensionalized and the most unstable mode was found numeri- 
meme by USe Of the initial value technique. The wind profile 
li a,U,sech* (y/1) was used. Simple CISK type heating was 
Mirnroduced tO détermine its effects on the growth of the 
waves. The jet wind profile was more unstable with easterly 
flow than westerly flow when no heating or friction was pres- 
meee nnen the jet profile was present only at the upper 
ievweie the ratio of the disturbance amplitude at the lower 
level to the disturbance amplitude at the upper level decreased 
toward the equator. When the heating rate was not large enough 
Memincrease the growth rate, the amplitude ratio remained small. 


mmcne heating rate becomes larger, the ratio increases to a 


fue larger than 1. 
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PrN EL RODUCT ION 


In the last tew years there has been growing evidence 
that barotropic instability may be a common occurrence at 
eitesZ200-mb level in the tropics. Krishnamurti [1970, 1971] 
Observed a peak at wave numbers 6-8 in the wind spectrum at 
this level which he proposed was caused by barotropic insta- 
bility. <A subsequent study by Kanamitsu, et al. [1972] showed 
Moae fOr the region between 15S and 15N waves of this scale 
memeecececive energy through barotropic interaction with the 
zonal mean flow and zonal wave number l. 

Mommie) examined the barotropic stability character- 
meeres Of 2 parallel] zonal current and found that instability 
mmposstole Only if the gradient of the absolute vorticity 
mmeoecs Sion somewhere in the region. The 200-mb level in 
mmemtrOpres is the outflow level of the giant monsoon circu- 
mreaon which is induced by differential heating (both sensible 
and latent) due to land-sea distribution. The basic flow 
memectrn Chus varies in the planetary scale. Detailed obser- 
vations have found that synoptic-scale moving disturbances 
memeeroelcvyet tend to form in the easterly jet south of the 
ec nigh {Krishnamurti, 1970] during the northern hemis- 
paere summer, This is a region of large vorticity gradient 
[iene tie mecessary condition for instability is occasionally 
mietstied locally. If these waves were the result of baro- 


meopic instability, they would extract energy from the mean 





flow and the long waves, since the long waves and the mean 
mow combine to give large vorticity gradients in these 
mee10ons. Ihe energy exchanges computed by Kanamitsu et al. 
[1972] are consistent with this view. 

Pits eeMestom ine barotropic instability of an upper 
ieee l easterly jet was examined. In addition to the calcu- 
ation of the growth rates, the downward propagation of energy 
was also investigated to examine the possibility that lower 
Mmeoposphere tropical disturbances may originate at the upper 
meres, olmple CISK type heating was introduced to ascertain 
how it affects the growth of the waves. 

Gimemstilay USCS a GuaSi-geostrophic two-level model. The 
equations are linearized and solved by the initial value 
Poomroue. Growth rate and wave structure are extracted after 


an exponential growth is established. 





ieee tii FORECAST EQUATIONS 


The quasi-geostrophic prediction equations were used as 
PiteetoOrecast Cuquations. It is recognized that the quasi- 
MeeetrOophic approximation 1S not very accurate near the 
Equator, but Matsuno [1966] has shown with a barotropic model 
that quasi-geostrophic motions are identifiable near the 
Pomator if the wavelength is not too long. Also, ina study 
with a baroclinic model, Yamasaki [1969] has found that the 
error is small at twenty degrees latitude. 

The quasi-geostrophic equations are simplified through 
Mmemise Of a two-level model. The basic model is very similar 
memene models used by Williams et al. [1971] and Robertson 
wee o|. the equations were linearized and solved numerically 
meeene anitial value technique. The growth rates and wave 
Pere Lures were then computed. 

ine Simole two-level model was constructed by dividing 
Bmemartmosphere into four layers with a constant pressure dif- 
memential of Ap/2 (Fig. 1), numbered 0 to 4 from top to bottom. 
The boundary condition at p = Pp, was wo = 0, which prohibits 
vertical propagation of energy. Charney [1969] has shown that 
such energy propagation is unlikely. At p = pe the boundary 
condition was w = We where Ww. was Pew eyret tonal ly induced a 
at the top of the Ekman layer. Charney and Eliassen [1949] 


derived the following expression for We 
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rb San 20 2 
a (2.1) 


_ Ae 
we = “°,8 (FR) 


where A, @oetne kinematic eddy viscosity, a is the surface 


mielow angle, and 2, is the surface geostrophic vorticity. 


Pp ewed 
2p ey —o L 
oe ee 
ree) 22> = Wa Peo e Se eose 2 
Bo, 5 Yi nnrcee---= 3 
a aT is aa 4 


Div iwo-lovye ll prediction model 


memias been shown by Holton et al. [1971] and Chang [1973] 
that this expression is reasonably valid as far equatorward 
Gemten degrees, but may be invalid close to the Equator south 


Gene JTCZ. The surface geostrophic vorticity was approximated 
by 
ty = 03 = V7 (2.2) 


iMiemavasi-pgeostrophic vorticity equation in pressure coordin- 


ates 1S 
OEE ; 2 EON pow 
moot KX VY + vCVv2¥) + oon earn CES) 


where ee = 2Qsin Do aem ZO Coys Wola, w = dp/dt, ¥ = stream 


O 
mime tion, is the reference latitude, and a is the radius 
Yo 


Paeene Carth. Observe that when w = 0, (2.3) becomes the 
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barotropic vorticity equation which is sufficient to describe 
mre barotropic instability. 


Popes equation at levels 1 and 3 gives 


r) . 2 Ona ee 
ac VY) OSE CY) 3 ae f. Ne ap (Zee 





ume * KXWY, + V(v7¥,) +8, SR - ee) - nee 
If it is desired to stop at some intermediate point, say the 
maepopause, rather than to include the entire atmosphere, 
pmemeter Layers could be used. 

leaweconsider sthie quasi-geostrophic first law of thermo- 


fmanmics in the form 


2 oly oly fo) _ 
dt dp IK XV¥ Urey a oe = Q (226) 
where 
Meela dl. ¢ 
oS me Las e oi ? (cz Ga) 


Seale is the temperature obtained from the standard atmosphere. 
The heating follows the formulation of Ooyama [1964], 

m@ermey and Eliassen [1964], and Kuo [1965], i.e., the con- 

densation heating is proportional to the convergence of moisture 


in the Ekman layer. The form used was given by Ogura [1964]: 


1 Cp Ts 30 
—2 nw, Ca) 


Peere ty 15 a non-dimensional parameter. This relationship was 


used in both rising and subsiding air representing a Fourier expansion 


LZ 





of the heating function. The factor 1/2 is therefore needed 
mommaintain the Correct amplitude of the heating perturbation. 


wicimenoimilcmappiecd atelevel 2 the result is 


5 ne APO2W. _ 
eee, - %3) * KXV(>-*)- vy, -¥,) - “PER = 0. (2.8) 


Define the following quantities 


¥, + ¥; 
uM = a 5 C29) 
Yn = 21th (2.10) 
which implies that 
Vo Ue 8 Tee C2320) 
ae (2 eZ 
ene Yep 1s proportional to the layer thickness and is therefore 


mrcasure Of the mean temperature in the layer. Using these 


Memmnitions, add (2.4) and (2.5) and divide the result by two, 





obtaining 
mon 42 : . 2 : 2 
ae My tT OK XV YY? WOVE) + (KX VY, + V(VEYD) 
aN fow 
M . LoWe _ 
1S cere Tee (CS) 


Memmeme second forecast equation, subtract (2.5) from (2.4) 


emeliminate w, by using (2.8) obtaining 
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Beep y? - 2 + |KX ; 2 2 2 
ee’ ur) Ya KX VY, * VCV* - We + IK X V¥,, ue) 
8 YT foWe 
Loa fap 1° (2.14) 
where 
ae 
2 a 2 1 4a 
N= po, ia 


miese are the prediction equations for the model. 


Pimeaough the normal mode method gives the phase speeds, growth 
miees, and wave structures for all modes, the initial value 
iapmeadch gives the phase speed, growth rate, and wave structure 
Seeonly the most unstable mode. This approach was chosen 
because only information about the most unstable mode was 
needed. 


teeseconvenient to write Vay and Ver in the following form: 


Vu Iemma. y,t) cos kx + B(y,t) sin kx, (Zo) 


Vr PG jemmey.,t) cos kx + D(y,t) sin kx. (2G) 


Gieme kK 1S the x wave number. The coefficients A through D 
feeeene Fourier amplitudes of the disturbance and E and F are 
Becemean zonal fields. 

iivese expressions for Vy and Wp are SMS tie ted into 
(2.13). The various sine and cosine terms are then separated, 
apematil products of the quantities A through D neglected. 


Equating the coefficients of the cosine terms gives 
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a ,d7A Seeoe Osea) 8 868°E. 68°F 
ot ay - AK) = Ki y dYay  oyoy? bys > ono 2 
JE OF 
K (sy B+ 57 PI] : pe IU ste). 
C2275) 
romechne Sine terms the result is 
ee Bon 2. = DEmonw or d-G . 0°5 a3F 
at (S52 BS is | dy oy? Cyd aya ° dy 3 . 
JE 0 F 
ket (AK ke te ire 
ae ) (B- D) 
(22-15) 


memeat the procedure for equation (2.14) for the cosine terms 


which gives 








Z 
we (Sy - CK? - cu?) = K[Ze 





dy? dy dy? y oy? dy 
era? = 2) Se oF Ry 
ay ay 3 ay 3 


Bo Dk as (K- 8) 3 - k*)(A-C). 


















(22419) 
Mmmemresults for the sine terms are 
p) ED oF 
3t (—> —Dk* - Dut) = k [5 LEC ae) oS (k? =) 
dE 37C ole 0 F 035 
By By2 ay ay2 * aysA* ays Ol 
fe B Ck + (K - s) (G5 -k?)(B- D), 
((27.2.0)) 
where 
_ fe gt VE 
K IRI Ge (2.204) 


IDS 





and 


Mi mecometuduhonal equation Lor vertical velocity 





4 n 96 Po 
: 8 re apes 
below: 
aT av 
ee te 
Apo ‘ot m9 
where 
noe 2 Ee 
mM fo oy 
oe oe 
Hl i) ox 





(——) 
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1S 


(2. 


(Ze 


(2. 


21) 


given 


22) 


23) 


24) 





IJI. BOUNDARY CONDITIONS AND FINITE 
re ee lNG oChHEME 


The finite difference scheme is illustrated below with 


a sample variable M: 


i 











M 
= or (54g 7 MG-4) > (B21) 
1° ae - 

a ae Magy - 2M, + My) (uel 
am. 1 _ : - 

mec epee ey GJ (UM, = OM; * Mi 2) I 


(55s) 
Peeeent 15 the grid index and H is the distance between grid 
points. 

Centered time differences were used for all quantities 
meee those involving friction and heating. The friction 
miomacating Cerms were computed at time (t- At). The first 
meen dail cases 1s a forward time step. The second order 
equations for time tendencies were solved by the exact method 
Seeeachtmyer {1957, p. 101]. 

mana and Nitta {1968a] studied the problem of using 
memes difference equations to solve dynamic instability 
peeiemems Of non-divergent barotropic currents with boundaries. 
In this study rigid boundaries are placed at y = + W/2 where 
wemiemenme tOtal width of the computational field. The boundary 


conditions were A= B=+Ce#D== Oaty = # W/Z. 


17 








OFZ 0.4 0.6 0.8 IO 


Fig. 2. Mean wind profile U. = a5Uy sech* (y/L) 
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IV. WIND PROFILE AND INITIAL CONDITIONS 


Wisestudy uses the following zonal wind profile: 


UL, = Bow seen Ey. 4.1 
: U, (y/L) (4.1) 


where 2 CimvceinOnreU dhe barotropic Stability of this 
Meorale has been analyzed by Lipps [1962]. Figure 2 contains 
mae profile (4.1) for as = wall 


The initial conditions are given below: 


A = sin ty/W (4.2) 
B = 0 (455) 
C = 0 aa) 
D= 0 (Aes. 
pee Fo Sr, ean. (y/L 4.6 
7) ia anh (y/L) (4.6) 

hea, = ay 
5 £oUGL tomes) 1) (4.7) 


Bemece tiiat the E and F fields are independent of time. 


19 





VV; NON-DIMENSIONALIZATION PROCEDURE 


ingmoncer to Simplify the interpretation of the results, 
eiemeduations are written in non-dimensional form as follows: 


foes, y = Ly , and t = L/t~. The stream function is 


Toma t ten 

ee Uy CO dy” * by (ey? t~) (5.1) 
and 

Ye = - LU sup (y)dy* + ¥p(x*,y7,t7) Ce) 


meme tnese expressions are substituted into (2.13) and 
(2.14), and when the products of the w's are dropped, the 


Hem@ecimensi1Oonal equations can be written: 











3 —o x * d2u OW 2 7 _ 
(se 3 Ux r)V* Way ¥ (B ae yx =e he Gy Ve) =) (S65) 
Pee ey G2 | » . St An 2 : 
(cers 2 ux) (V E) pep + (8 dy?) x Vets or) Vv (Way wer 
(34743 
where 
rn -1 ell Sj 2 A 2 
= 22 as 2 A pe eas: e 
€ l= Vo ? B Bot U ; r = 2Ap/P.i£ (7F) (So) 
and 
Chee : 
Aan nmax L pogt’ |pplsin 2 o as ‘ (5.6) 
Uo4t fo J, at | 


Mmiemouantity € 1S a rotational Froude number, 8 is a non- 
dimensional form of 8, Tf is the ratio of the square root of 


the Ekman number [Greenspan, 1968] to the Rossby number, and 
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6 is the heating parameter. The primes on the independent 
memtables have been dropped. These equations are equivalent 


Meme set 2.17, 2.18, 2.19, and 2.20 which were solved 


numerically. 


2) 





Vie COMPULATTONAL PROCEDURE 


The following constants were used in the computations: 


De = 100 centibars, 

hs = igi ometme ies, 

W ae SOLOL0 EE et ane) Kono oe ee 

Uo —wUeete rom ic mes ceOnd 

LS = varied with latitude, 

O So UsMecrems per Second Squared per 


centibar squared, 


B Seucercan nthe lata tude , 

Ap = 50 cCentibars, 

Ae = 10 meters per second squared, 
OL ee edconces . 

Wie = varied with latitude 

H = 50 kilometers, 

ie Sele OU OUTS. 

“max uO 


d6/ap= -0.6 degrees per centibar. 
A forecast period of 100 days was used so that the most unstable 
mode would have reached the maximum growth rate. For conven- 
ience, growth rates were computed over the last two days of 
[miemrorecast period. This was done by assuming the amplitudes 
Gould be written as 

econ U/L) t (cea) 


x) = @. S90 (WB a (6.2) 


ee 





where n is the non-dimensional growth rate and ao 1Sa CONS tane, 


By forming the ratio x,/x,, the growth rate can be shown to be 


Seen (x5 / xX, ) 


Zo 





Vilieeke oul ls 


A series of experiments were conducted over a wide variety 
Of conditions. The more important experiments are given in 


ites! as a function of a,, a,, B, ce, 6, T, n, and k*, where 


3? 
Mmereoresents the maximum growth rate for each set of conditions 
and k~ is. the corresponding non-dimensional wave number kL. 
When some of these parameters were computed over a range of 
values, the range is given. 

Experiment 1 was conducted to determine the difference in 
growth rate characteristics between easterly and westerly jets. 
iments experiment no heating or friction was included and the 
mean flow was the same at both levels. In order to obtain a 
Continuous variation between the easterly and westerly flows, 
the non-dimensional parameter § (see equation 5.5) was defined 
to assume the sign of U. Thus 8 is positive for westerly flow 
omeemegative for easterly flow. In the remaining experiments 
the jet will be easterly, but 8 will be regarded as positive. 
For each 8 the maximum growth rate and the corresponding wave 
Meer were determined. Figure 3 contains these quantities 
as a function of 8; the maximum growth rate of n = 0.1888 and 
occurs at 8B = -0.352. The corresponding wave number is k* = 
0.99. This value of k” is close to the minimum k° which occurs 
when B = -0.112. Poems hows cleanmby that the easterly jet 
aS more unstable than the westerly jet. In particular, the 


Mmuowth rates for the easterlies over a range of beta are 
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Seceateretnan the value for beta = 0. This is in contrast to 
Biemenormmlal notion that the presence of beta has a stabilizing 
emerect fOr inStabilities. 

Experiment 2 investigated the downward propagation of 
Piereyetrom a barotropically unstable disturbance as a func- 
tion of the latitudinal location of the Wee, AM Sele eke 
faerntroduced at the upper level and the lower level mean 
meee was set to zero. Surface friction was included, but 
heating was neglected. <A measure of the downward propagation 
/|¥, | 


fmeemmcd Lrom the exponentially unstable solution. Figure 4 


is the ratio R = |¥,| where these amplitudes are 


maxX max 


Contains isolines of R as a function of k”~ and latitude. The 

curve no, Shows the wave number which has the maximum growth 

meow at Cach latitude. The R value decreases as one moves 

inom high to low latitude as predicted by Charney [1963]. 

The R value also decreases with increasing wave number at 

maamebatatude. This figure shows the dependence on epsilon. 
Poperimemes 3, 4, and $§ all contain heating and were con- 

ducted with a fixed latitude of 10 degrees, otherwise they 

are the same as experiment 2. In these experiments a weighting 

factor was applied to the heating following Bates [1970] and 


Yamasaki and Wada [1972]. The following expression was used 


- 2 
n Taye COM yi We (feels) 


This form implies that moisture is available for the heating 
in the region |y|<L. If this was not done, CISK heating would 


Seeumein regions away from the basic jet. In experiment 3 
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where ee 4 the non-dimensional heating parameter is 

pee 0.165.) Freure 5 shows the growth rate (n) and the ampli- 
tude ratio (R) as a function of the non-dimensional wave number 
foe «=Note that the maximum growth rate of n = 0.187 occurs at 
k~ = 1.0. These two values are very close to the values ob- 
mamed Without Weating. The R is fairly large at the lower 
values of k*“ with a maximum at 0.5. R becomes very small in 
the area of maximum growth. Apparently the barotropic insta- 
bility effects dominate in the areas of maximum growth while 
mfesieating 1s relatively more important for k”~ < 0.5 as 


eeeraenced by the larger values of R. This effect is even 


mere pronounced in experiment 4. 


iy cxpelmmunemuns where es Gwe Mom aimenstonal heat ine 
parameter is 6 = 0.247. Figure 6 shows that the maximum growth 
meee OL Th — 05186 occurs at kX = 1.0 just as in experiment 3. 


As k”~ varies away from 1.0 to either extreme the growth rate 
meeoaches 0el0/. The R is larger in both of these regions 
mimere the largest value occurs for the largest value of k”. 
Peeoccoldary maximum in R occurs at k~ = 0.6. This indicates 
that heating dominates the solution in the regions which have 
mie Larger values of R. 


mone eniment 5S where np 8 the non-dimensional heating 


max 
Parameter is §6 = 0.330. Figure 7 shows that the growth rate 
ieedarly constant over the k~ considered. This is similar 

to the results obtained by Williams and Robertson [1973] who 


acquired essentially a constant growth rate with a hyperbolic 


tangent wind profile using larger values of 64. 
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The growth rate n (solid line) and the corre- 
sponding wave number k’~ (dashed line) as functions 
eee Ommeperiment sl. Westerly rlow occurs for 
peu edndmeasterly flow for 6 < 0. 
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Experiment 3. 
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Vor CONCEUSIONS 


Recent studies suggest that disturbances may form through 
mweceropic Instability in the easterly jet south of the 
i@eetan High at 200-mb. In this thesis the instability of an 
ematytical jet profile was investigated with a two-level 
quasi-geostrophic model. The downward propagation of energy 
mean dad barotropically unstable upper level disturbance was 
merestigated to see 1f it could initiate lower tropospheric 
Ges turbances. 

The hyperbolic secant squared wind profile was used in 
all of the numerical experiments. The first experiment, which 
@omeained no heating or friction, compared the growth rates 
With easterly and westerly flow. The easterly jet was the 
more unstable and its maximum growth rate was well above the 
growth rate for ~ = 0. The energy propagation was studied 
by introducing the jet at the upper level with no mean flow 
meow. ihe ratio (R) of the disturbance amplitude at the 
[imeem ievel to the disturbance amplitude at the upper level 
Peemecompuued. With surface friction included, R decreased 
from high latitudes to low latitudes; at each latitude R 
decreased with increasing wave number. When CISK type heating 
Poemeeneroduced With a large enough heating coefficient, the 
Puewemerate became independent of wave number. In this case 
R had a finite value and it increased with increasing wave 


Moet. For smaller values of the heating coefficient the 


J 





Peorhonutcmmlseanpriity Controls the growth rate and R remains 
Merny Small. 

The result that the easterly jet 1S more unstable than 
the westerly jet supports the view that barotropic instabil- 
[ir nOccurs souti or tne fibetan High. The small value of the 
Mmplrtide fatio at low latitudes agrees with Charney's [1963] 
conclusion that the vertical energy propagation is small. 

The inclusion of CISK heating did increase the low level 

Mime rTtudce,spuemtmrs Could not occur unless a finite amplitude 
oe Ppresehumdeethie Lower level. Since these results were 
Meeaiiied Witiea GUuasil-peostrophic model with only two levels, 
if iS important to carry out further studies using more 
levels and also allowing for gravity waves. This would 
meouwire a multi-level primitive equation model in either 
Spger eal coordinates or on an equatorial 6-plane. <A more 
sophisticated heating parameterization, such as that by 


Arakawa and Schubert [1974], should be used. 
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